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Review 

From computational modelling of the intrinsic 
apoptosis pathway to a systems-based analysis of 
chemotherapy resistance: achievements, perspectives 
and challenges in systems medicine 

ML Wurstle^ E Zink^ ^ JHM Prehn^'^ and M Rehm*'^'^ 

Our understanding of the mitochondrial or intrinsic apoptosis pathway and its role in chemotherapy resistance has increased 
significantly in recent years by a combination of experimental studies and mathematical modelling. This combined approach 
enhanced the quantitative and kinetic understanding of apoptosis signal transduction, but also provided new insights that systems- 
emanating functions (i.e., functions that cannot be attributed to individual network components but that are instead established by 
multi-component interplay) are crucial determinants of cell fate decisions. Among these features are molecular thresholds, 
cooperative protein functions, feedback loops and functional redundancies that provide systems robustness, and signalling 
topologies that allow ultrasensitivity or switch-like responses. The successful development of kinetic systems models that 
recapitulate biological signal transduction observed in living cells have now led to the first translational studies, which have 
exploited and validated such models in a clinical context. Bottom-up strategies that use pathway models in combination with 
higher-level modelling at the tissue, organ and whole body-level therefore carry great potential to eventually deliver a new 
generation of systems-based diagnostic tools that may contribute to the development of personalised and predictive medicine 
approaches. Here we review major achievements in the systems biology of intrinsic apoptosis signalling, discuss challenges for 
further model development, perspectives for higher-level integration of apoptosis models and finally discuss requirements for the 
development of systems medical solutions in the coming years. 
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Facts 

• Mathematical modelling has provided mechariistic insight 
into cell death decisions during intrinsic apoptosis signal- 
ling by identifying network-level control functions that 
emanate from multi-protein interplay. 

• Validated systems models of intrinsic apoptosis can 
accurately predict signal transduction in individual cell 
lines, kinetically and quantitatively. 

• In translational studies, systems models of intrinsic 
apoptosis show potential for predicting tumour responsive- 
ness to genotoxic chemotherapy. 

Open Questions 

• Can systems models that accurately predict the respon- 
siveness of cell lines to targeted apoptosis inducers and 
sensitisers be validated clinically? 



• Will we be able to develop systems-based diagnostic 
tools that in the future may contribute to the development 
of personalised medicine and innovative patient 
stratification? 

• Can clinical research environments adapt to fulfil the needs 
of translational systems biology and systems medical 
research? 

The intrinsic apoptosis pathway can be initiated by a multitude 
of diverse stress signals (Figure la), including DNA damage, 
ER stress, cytoskeletal destabilisation, kinase and protea- 
some inhibition. These stresses result in the transcriptional or 
post-translational activation or accumulation of BH3-only 
proteins, pro-apoptotic members of the B-cell lymphoma-2 
(Bcl-2) protein family. BH3-only proteins activate the 
proapoptotic proteins Bax and Bak and/or antagonise 
anti-apoptotic members such as Bcl-2, Mcl-1, Bcl-xL. 
Activated Bax and Bak can homo-oligomerise to form pores 
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Figure 1 Convergence of complex apoptosis signals into a cell death decision - a simplified overview, (a) The upstream sensing of cellular stress and damage and the 
resulting initiation of apoptosis signalling is subject to a complex interplay of multiple signalling pathways, (b) Stress signals are integrated at the level of Bcl-2 family proteins. 
Their interactions and the regulation of their interactions dictate whether cells can undergo MOMP. (c) The apoptosis execution phase entails a number of positive feedback 
loops which result in the amplification of effector caspase activities and thereby ensures efficient cell death execution. Negative regulations in (b) and (c) were not visualised to 
increase clarity 



in the outer mitochondrial membrane. Through these pores, 
the pro-apoptotic proteins cytochrome-c and Smac are 
released into the cytosol. The complex interactions of Bcl-2 
family protein members therefore can integrate highly diverse 
upstream stress signals into a response of mitochondrial outer 
membrane permeabilisation (MOMP) (Figure 1b). In the 
cytosol, cytochrome-c triggers apoptosome formation and 
caspase-9 activation, whereas Smac antagonises caspase- 
inhibiting lAP proteins such as XIAP. Activated caspase-9 
proteolytically activates effector caspases-3 and -7. These 
effector caspases initiate feedback amplification loops and 
are the ultimate executioners of apoptotic cell death 
(Figure 1c). The execution phase is a largely dormant 
sub-system of apoptosis signal transduction, whereas the 
signalling network upstream of BH3-only protein signalling 
exists in a dynamic steady-state and continuously responds 
and adapts to stress signals. 

Owing to the convergence of diverse stress signalling 
pathways into a common execution phase, and the fact that 
biochemical data of proteins and reactions involved in this 
phase have been made available, the execution phase has 
been the first target of systems biological studies that aimed to 
identify features, which determine and control cell death 
decisions. More recently, the first detailed systems 
approaches towards the understanding of Bcl-2 protein 
interactions and MOMP have been developed. In contrast, 
the very early processes that initiate the intrinsic apoptosis 
pathway, such as the transcriptional control of BH3-only 
protein production or their post-translational activation, have 
remained rather elusive or inaccessible to systems biological 
studies so far. An exception to this is the extrinsic apoptosis 
pathway that requires the post-translational activation of the 
BH3-only protein Bid (reviewed by Lavrik and colleagues in 
this special issue). 

In the following, we will briefly describe the challenges faced 
for the systems biological analysis of early apoptosis initiation 
events and will then discuss and critically evaluate systems 
biological studies focusing on Bcl-2 family protein interactions, 
MOMP regulation and apoptosis execution. Recent reviews 



on mathematical modelling techniques used in systems 
biological studies of apoptosis signalling may serve to obtain 
deeper methodological insight in this context. ^"^ Open 
questions and important aspects to be considered for future 
systems-level investigations of apoptosis signalling will be 
covered in this context as well. We review studies that 
successfully applied apoptosis systems models in transla- 
tional and clinical research settings, and will provide perspec- 
tives for the higher-level integration of apoptosis pathway 
modelling to allow multi-scale modelling beyond the 
single-cell level. Finally, we will highlight requirements and 
challenges faced for the clinical evaluation and implementa- 
tion of apoptosis systems biology within the context of the 
newly emerging field of systems medicine. 

The Initiation of tlie Intrinsic Apoptosis Pathway - 
A Black Box for Kinetic Systems Modelling 

The initiation of the intrinsic pathway currently remains largely 
uncharted territory for kinetic systems modelling. This is 
partially owed to the diversity of stress signals that can initiate 
intrinsic apoptosis, the large number of BH3-only proteins, 
their possible combinations and their relative contributions 
towards signal transduction, and the respective mechanisms 
controlling the production or activation of individual BH3-only 
proteins (transcription, translation, post-translational modifi- 
cation). In addition, experimental tools that allow the collection 
of quantitative kinetic data for these processes at the single- 
cell level and which are required both for model parameter- 
isation and for experimental validation of model predictions 
have not yet been developed. However, gene transcription 
and protein expression dynamics of other cellular signalling 
processes have been successfully monitored in real time in 
recent years."^"^ Such tools could therefore be adapted to 
investigate early stages of apoptosis initiation at the single-cell 
level. The need for a single-cell-based analysis is of particular 
importance for the initiation phase, as the majority of cell- 
to-cell heterogeneities arise from variability in gene transcription 
and protein expression that affect early apoptosis signalling.^"^ 
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Besides expression noise in apoptosis signalling proteins, 
interplays with non-apoptotic signalling pathways may also 
contribute significantly during early stages of signal transduc- 
tion. For example, proliferation signalling is known to negatively 
regulate apoptosis susceptibility, and pathways involved in 
sensing and counteracting the damage inflicted by inducers of 
intrinsic apoptosis (e.g., genotoxic agents and ionising irradia- 
tion) contribute towards defining resistance towards apoptosis 
initiation. Large-scale statistical systems modelling building 
on complex data sets generated by high content and high- 
throughput experimental approaches rather than detailed 
kinetic pathway models therefore currently seem to be 
superior in describing and tackling the complexity of these 
pathway crosstalks. Indeed, such studies have successfully 
identified network components, which regulate and contribute 
towards defining apoptosis responsiveness.''^"'"^ 

Systems Analysis of the Regulation of IVIOMP 

During intrinsic apoptosis, MOMP is essential for triggering 
the apoptosis execution phase. MOMP is regulated by the 
family of Bcl-2 proteins. So far, more than a dozen members of 
this protein family have been identified, and these engage in a 
complex interplay of activation, inactivation and multimerisa- 
tion.''^ All Bcl-2 family members contain at least one of four 
different Bcl-2 homology (BH) domains (BH1, BH2, BH3 and 
BH4), which define their functional role as either being pro- or 
anti-apoptotic. The anti-apoptotic Bcl-2 family members 
(Bcl-2, Mcl-1, Bcl-xL, A1 and Bcl-w) contain all four BH 
domains. They prevent apoptosis induction by antagonising 
the pro-apoptotic Bcl-2 proteins. The pro-apoptotic fraction of 
the Bcl-2 family can be divided into BH multi-domain proteins 
(Bax, Bak) and BH3-only proteins. All BH3-only proteins can 
antagonise anti-apoptotic Bcl-2 proteins, whereas the BH3- 
only proteins Bim, Puma and Bid can also directly activate Bax 
and Bak (Figure 2a). Individual BH3-only proteins or 
combinations thereof are induced or activated in response 
to stress signals that initiate the intrinsic apoptosis pathway, 
such as irradiation, genotoxic drugs or endoplasmic reticulum 
stress. Through Bax/Bak pores, pro-apoptotic factors are 
released into the cytosol to trigger apoptosis execution. 
Therefore, the interplay of Bcl-2 family proteins that controls 
MOMP is highly complex. 

Mathematical modelling studies of varying levels of 
complexity and detail have provided important insight into 
systems-emanating features that arise from Bcl-2 proteins and 
have helped to mechanistically explain experimentally 
observed signalling characteristics. This includes the rapid 
switch-like kinetics of MOMP, the robustness of the ultra- 
sensitive pore formation process to variability in protein 
expression levels, as well as the spatiotemporal spread 
of MOMP signals because of reaction-diffusion signalling. 
Such studies made use of kinetic models based on ordinary 
differential equations''^ ''^ or spatial modelling approaches 
such as stochastic cellular automata simulations and partial 
differential equations based models. ^'^^'^^ Following from 
initial systems theoretical investigations that studied the 
systems-level consequences, which arise from the known 
interaction topologies, these mathematical modelling 
approaches subsequently were successfully linked with 
quantitative experimental data. 



Experimental observations of MOMP at the single-cell level 
demonstrated that the release of mitochondrial inter- 
membrane space proteins is a rapid process that is completed 
within less than 5min.^° Only when investigated at a 
resolution of seconds, stimulus-specific differences in MOMP 
kinetics become apparent.'' Given that the overall duration of 
apoptosis signal transduction from initiation to execution 
typically requires many hours, the process of MOMP is a 
remarkably rapid unidirectional switch-like response of 
membrane permeabilisation. One of the most central systems 
biological problems therefore is to understand how this 
permeabilisation switch can emanate from the interplay of 
Bcl-2 family proteins. In this context, experimentalists have 
discussed at length whether enabler and activator BH3-only 
proteins can drive rapid Bax pore formation equally efficiently 
(indirect versus direct activation scenario). In the indirect 
activation scenario, anti-apoptotic Bcl-2 family proteins 
antagonise active Bax or Bak proteins. Enabler BH3-only 
proteins interfere with this interaction and liberate active Bax 
and Bak by binding to anti-apoptotic Bcl-2 family members. In 
the direct activation scenario, the BH3-only protein members 
are divided into enabler and activator proteins. Although 
enabler proteins (such as Bad, Bmf, Noxa and Bik) only target 
anti-apoptotic Bcl-2 family proteins, activator proteins (such 
as Bim, Puma and Bid) interact with and directly activate Bax 
and Bak, thereby triggering MOMP. Several systems theore- 
tical studies investigated the plausibility of these two activa- 
tion models. Chen et al?^ tackled this question long before 
conclusive experimental evidence was available, using a 
continuous deterministic modelling approach based on 
ordinary differential equations. A rigorously simplified interac- 
tion network considered only one species each of a Bax-like, a 
Bcl-2-like and an enabler and activator BH3-only protein. 
Parameter screening showed that ultrasensitivity in pore 
formation (i.e., the rapid formation of large amounts of Bax 
pores upon reaching a critical threshold of BH3-only proteins) 
is very unlikely to manifest from an indirect activation model. 
In contrast, a model of the direct activation scenario displayed 
ultrasensitive Bax pore formation processes for a large 
number of parameter sets. On the basis of this study, Bax 
pores would be expected to form in a switch-like manner once 
a critical input of activator BH3-only protein has accumulated. 
In contrast to the indirect activation model, the direct activation 
scenario furthermore is robust to perturbations in the 
respective protein amounts. These findings were largely 
validated in a recent study, which provided the so far most 
detailed mathematical model of the Bcl-2 interaction network, 
including a biologically meaningful quantitative parameterisa- 
tion of protein amounts and rate constants^'' (Figures 2b and c). 
In particular, it was shown that an indirect Bax/Bak 
activation mechanisms through enabler BH3-only proteins 
require strong antagonising regulatory mechanisms such as 
very high Bax/Bak pore degradation kinetics to prevent basal 
rates of spontaneous pore formation.^'' In contrast, a direct 
activation scenario results in ultrasensitive pore formation 
responses that are robust to cell-to-cell heterogeneities in 
protein expression amounts and that were in line with 
previously published experimental data^'' (Figure 2b). 
Integrating both activator and enabler BH3-only proteins in 
this model was sufficient to provide very good predictions on 
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the MOMP susceptibility and cell death responsiveness of a 
number of cancer cell lines treated with genotoxic drugs'*^ 
(Figure 2c). The kinetics of Bax pore formation could also be 



reproduced qualitatively using a stochastic cellular automaton 
modelling, which represented Bax pore formation in 
membrane patches.''^ This methodological approach was 
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subsequently combined with quantitative experimental single- 
cell measurements in which absolute cytosolic and mito- 
chondrial Bax quantities were measured in real time during 
MOMP.^^ Not only could the mathematical model reproduce 
pore formation events at the experimentally measured 
amounts of mitochondrially accumulated Bax (Figure 2d), it 
was also found that only a small fraction of the entire Bax pool 
(-^5%) needs to be activated in presence of activator and 
enabler BH3-only proteins to trigger pore formation. 
Mechanistically, it appears that the unidirectional switch-like 
process of MOMP emerges as a consequence of cooperative 
Bax oligomerisation, which occurs once a critical amount of 
active Bax has accumulated.^^ MOMP then rapidly spreads 
through the cell body, a process that kinetically can be 
re-modelled by considering BH3-only proteins distributing 
spatiotemporally by reaction-diffusion processes.^ Interest- 
ingly, the velocity of these MOMP waves may depend on the 
rate of BH3-only protein production. 

In agreement with the mathematical simulations, the 
importance of the direct activation mechanism for rapid pore 
formation has been supported experimentally by both 
mechanistic studies and by measurements of relative 
protein abundances of pro- and anti-apoptotic Bcl-2 family 
proteins. ^^'^^'^^ Nevertheless, to fully resemble the biological 
characteristics of MOMP, contributions from enabler BH3- 
only proteins towards Bax/Bak pore formation need to be 
taken into account as well.^ '''^^'^^ Although pore formation is 
an ultrasensitive, non-linear response, it needs to be noted 
that also the prevention of pore formation by anti-apoptotic 
Bcl-2 proteins is controlled in a highly non-linear fashion. 
Gaudet et al.^^ showed by both mathematical modelling and 
experimental measurements that critical thresholds of Bcl-2 
exist beyond which cells suddenly lose their capacity to 
undergo MOMP (Figure 2e). Importantly, this transition occurs 
across a narrow protein concentration range, and small cell- 
to-cell differences in Bcl-2 protein expression therefore can 
contribute towards explaining that only fractions of HeLa cell 
populations may respond to apoptotic signals.^"^ 



Open Questions and Considerations 

The mechanistic and quantitative description of the individual 
interactions controlling Bax/Bak pore formation and MOMP 
kinetics is still far from complete. Especially, the complexity of 
the multiple interactions between the subfamilies of the Bcl-2 
proteins poses a challenge. Also, interactions appear to take 
place with significantly different reaction rates in aqueous and 



membranous environments, and contributions from non-Bcl-2 
proteins may need to be considered as well.''^'^^"^^ Additional 
studies aimed at determining binding characteristics of Bcl-2 
proteins at or in membranes are therefore required, and 
systems biological studies on the control of MOMP are likely to 
significantly benefit from such data and will allow a deeper 
biological understanding of MOMP regulation. Furthermore, 
current models have not taken into account the continuous 
cycling of Bcl-2 proteins between the soluble cytosolic fraction 
and membrane-associated pools of the proteins, a process 
recently described for the Bax/Bcl-xL interplay.^^'^° In addi- 
tion, various Bcl-2 proteins were shown to exhibit functions 
beyond apoptosis induction or prevention. For example, anti- 
apoptotic Bcl-2 proteins suppress autophagy by binding and 
inhibiting Beclin-1, a key protein required for early steps of 
autophagosome formation.^'' Bcl-2 proteins also co-regulate 
mitochondrial fusion and fission as well as mitochondrial 
bioenergetics.^^'^^ Apoptosis systems modelling could there- 
fore benefit from integrating such protein multi-functionalities 
as well as existing metabolic models in order to determine 
apoptosis susceptibility also in the context of cellular energy 
supply and mitochondrial fitness. 



Systems Analysis of tlie Apoptosis Execution Pliase 

As a consequence of MOMP, cytochrome-c and Smac are 
released into the cytosol. Cytochrome-c, together with dATP/ 
ATP, binds to monomeric Apaf-1. Apaf-1 then undergoes a 
conformational change and oligomerises into a stable 
heptamer that recruits procaspase-9. Procaspase-9 is acti- 
vated upon recruitment and proteolytically processed. The 
main procaspase-9 substrates are the zymogens procas- 
pase-3 and -7. Caspases-9, -3 and -7 can be inhibited by 
XIAP, which itself is inhibited by Smac. Active caspase-3 is a 
key driver for rapid apoptosis execution by cleaving hundreds 
of cellular substrates and directly or indirectly causes the 
characteristic morphological changes associated with apop- 
totic cell death, such as cellular shrinkage, membrane 
blobbing and DNA f ragmentation^^'^^ (Figure 3a). The kinetics 
of apoptosis execution resulting from such a linear sequence 
of events would be expected to be strongly sensitive to 
changes in the relative abundance of the involved proteins. 
However, real-time measurements in living cells demon- 
strated the opposite in that apoptosis execution in various cell 
lines with demonstrably different protein expression profiles 
proceed with nearly identical kinetics. ^'^'''^^ Furthermore, 
substrate cleavage by effector caspases resembles an 
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Figure 2 Systems analysis of Bcl-2 protein signalling and MOMP. (a) Simplified overview of Bcl-2 protein interactions. Associated systems-level features 
and underlying molecular causes are listed, (b) Mathematical simulations of the calculated BH3-only protein amounts required to induce MOMP for scenarios of direct and 
indirect Bax/Bak activation. A number of cell lines known to be capable of undergoing MOMP would require excessive amounts of BH3-only proteins in an indirect 
activation setting. For the direct activation scenario, biologically justifiable amounts of BH3-only protein inputs were calculated for all cell lines, demonstrating that the direct 
activation model is in line with experimental evidence. Calculations were performed using kinetic models of Bcl-2 protein interactions. ""^'^^ (c) Amounts of BH3-only 
proteins required for MOMP in a direct activation model were calculated for a number of colorectal cancer cell lines. Calculated BH3-only amounts correlate well with the 
responsiveness of the cell lines to genotoxic apoptosis inducers. Calculations were performed using a kinetic model of Bcl-2 protein interactions.''^ (d) Stochastic 
simulations of Bcl-2 protein interactions and Bax pore formation in membrane patches.^^ Modelling results indicate that the initiation of Bax pore formation correlates well 
with experimentally measured amounts of activated Bax at mitochondria, (e) Simulations of the effect of Bcl-2 protein expression amounts on the cellular competence to 
undergo MOMP resemble experimentally observed cellular behaviours. With increasing Bcl-2 protein amounts, cells undergo a sharp transition from MOMP competent to 
MOMP resistant. Simulations were conducted using a kinetic model of apoptosis signal transduction.^"^ Display items in panels (b-e) were reproduced or modified from 
previous publications^^'^^'^^'^^ 
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all-or-none process, which is independent of the stimulus and 
stimulus strength and which is triggered by rapid effector 
caspase activation.^ The robustness of the substrate 



cleavage switch is established by amplifying feedback loops 
within the apoptosis execution network.^^ The main positive 
feedbacks are the caspase-3-mediated removal of the 
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lAP-binding domain of caspase-9, which efficiently alleviates 
the inhibition of caspase-9 by XI AP,'^^ and the sequestration of 
XIAP away from active caspase-9 by active caspase-3.'*'' The 
latter positive feedback manifests prominently when the 
amount of caspase-3 exceeds that of caspase-9 and XIAP, 
as at these conditions activated caspase-3 can serve as a sink 
that depletes XIAP and thereby re-activates caspase-9. It was 
also shown that caspase-3 may process and activate its own 
pro-form, albeit with significantly lower catalytic activity.^^ 
A detailed mathematical model based on ordinary differential 
equations was developed to investigate the execution signal- 
ling network in HeLa cervical cancer cells. Parameterised 
with experimentally determined protein amounts, the mathe- 
matical simulations matched the experimentally recorded 
kinetics of effector caspase-dependent substrate cleavage 
very accurately (Figure 3b). Simulations were also sufficiently 
accurate to predict the time required to activate effector 
caspases in various other widely used cell line systems, 
including colorectal and breast cancer cell lines^^ (Figure 3c). 
The in silico analysis of the consequences of altering the 
amounts of procaspase-9, Smac, Apaf-1 and procaspase-3 in 
HeLa cells indicated that the execution network remains 
functional and assured effector caspase activation across a 
wide range of protein concentrations. However, further 
analysis showed that the elevation of XIAP above a critical 
threshold concentration would result in a sudden collapse of 
the all-or-none response and at such conditions promote a 
limited and sublethal activation of effector caspases^^ 
(Figure 3d). Importantly, the transition between efficient and 
sublethal induction of the apoptosis execution phase occurs 
across a narrow range of XIAP concentrations, indicating that 
modulations of XIAP amounts can transfer the network from 
an apoptosis competent to an apoptosis resistant state 
(Figure 3e). A systems theoretical investigation of the 
signalling network also identified that the execution network 
presents bi-stability characteristics in which the system can 
undergo a transition from a 'live' state to an irreversible state of 
lethal effector caspase activation.'^'' It is important to note that 
the switch-like transition from apoptosis competent to 
resistant is not a function of XIAP alone but rather emanates 
from the interplay of multiple network components. This 
transition is physiologically relevant and has been observed 
during the differentiation of post mitotic cells. For example, 
neurons and cardiomyocytes become apoptosis resistant by 
re-adjusting the relative abundance of execution network 
proteins during differentiation. "^^""^^ Subsequent to quantita- 
tive kinetic simulations of apoptosis execution during intrinsic 
apoptosis, these signalling processes were also integrated 
successfully into models of extrinsic MOMP-dependent 
apoptosis signalling."^^ 



Open Questions and Considerations 

Although experimentally validated models of apoptosis 
execution have been remarkably accurate for the investigated 
cell line systems, they are of limited detail concerning the early 
molecular processes following MOMP. Especially, the forma- 
tion of the Apaf-1 *caspase-9 apoptosome and caspase-9 
activation so far have been implemented mostly as input 
functions that remodel the experimentally determined in vitro 
kinetics of apoptosome formation and caspase-9 activation. 
Although separate theoretical investigations into apoptosome 
formation were conducted, these only superficially integrated 
the available information on apoptosome formation mechan- 
isms."^^'"^^ A more detailed systems-level analysis would also 
be required to investigate and understand the potency of 
processes that were recently described to suppress Apaf-1 
heptamerisation and caspase-9 activation. For example, 
phosphorylation of Apaf-1 by 90-kDa ribosomal S6 kinase 
(RSK) prevents Apaf-1 oligomerisation."^^ RSK is activated by 
the mitogen-activated protein kinase pathway, which is 
frequently constitutively (over-)activated in cancer cells. ^° 
If Apaf-1 is phosphorylated by RSK, it recruits the adaptor 
protein 14-3-3s, and 14-3-3s then blocks the binding of 
cytochrome-c, Apaf-1 activation and apoptosome forma- 
tion."^^ Similarly, Cdc6, a AAA+ ATPase that is upregulated 
in various cancers and linked to carcinogenesis,^'' prevents 
apoptosome formation by binding to cytochrome-c-activated 
Apaf-1 .^"^ Of note, caspase-3 cleaves and inactivates Gdc6 
during apoptosis execution. This feedback could free Apaf-1 
from Gdc6-dependent inhibition and may restore apoptosome 
formation and thereby could enhance effector caspase 
activation. An extension of existing mathematical models 
could therefore serve to identify the context-dependent 
potencies of these mechanisms in preventing apoptosis, 
and the likelihood and robustness of inefficient and sub-lethal 
caspase activation. Current models also do not take into 
account that the apoptosome acts as a proteolytic-based 
molecular timer.^^'^^ Once procaspase-9 is activated at the 
apoptosome, it undergoes auto-processing. Owing to its weak 
association with the apoptosome backbone, processed 
caspase-9 is released into the cytosol and thereby becomes 
inactivated. At which conditions this molecular timer substan- 
tially contributes to cell death decisions has not been studied 
so far. 

Owing to multi-protein reaction-diffusion processes during 
apoptosis execution, spatial anisotropies observed during 
the upstream event of MOMP are efficiently eliminated 
during apoptosis execution in epithelial cells, as was shown 
by theoretical and experimental analysis. However, 
spatially coordinated and limited apoptosis execution has 
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Figure 3 Systems analysis of the apoptosis execution phase, (a) Simplified overview of core signalling processes during the apoptosis execution phase. Associated 
systems-level features and underlying molecular causes are listed, (b) Kinetic simulations of substrate cleavage by effector caspases correlate well with experimental 
measurements in individual parental HeLa cells.^^ (c) Mathematical calculations of the time required from IVIOIVIP until effector caspases activation in different cancer cell lines 
correlate well with experimental measurements.^^ (d) Predictions on the consequence of XIAP overexpression on substrate cleavage by effector caspases were 
experimentally confirmed by measurements in HeLa cells overexpressing XIAP. Substrate cleavage is prolonged, submaximal and sublethal.^^ (e) Mathematical analysis of 
the consequences of increasing XIAP amounts on effector caspase activity in HeLa cells. The simulations highlight that cells can undergo a sudden transition from apoptosis 
competent to apoptosis resistant once a critical XIAP concentration is exceeded. Display items in panels (b-e) are results from kinetic systems models of apoptosis execution 
and were reproduced or modified from previous publications^^'^^ 
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been observed in neuronal cells. ' A spatiotemporal 
systems analysis for such large, outstretched cells 
may therefore be required to understand whether our 
current knowledge on the reaction-diffusion properties of 
apoptosis execution proteins can explain the experimentally 
observed features of localised activation of effector 
caspases. 



Translational Application of Systems IVIodels of tlie 
Intrinsic Apoptosis Pathway 

Although defects in apoptosis signalling are a hallmark of 
cancer,^"* associating the expression of individual apoptosis 
proteins with cancer therapy resistance so far largely failed. 
It is therefore hoped that pathway systems models, which 
take into account multi-protein interactions, may be superior 
predictors for the responsiveness of cancer tissues to 
apoptosis-inducing treatments. First proof-of-concept studies 
that evaluated the translational capacity of apoptosis systems 
models have been conducted in recent years. Mathematical 
modelling of the Bcl-2 family network, parameterised with 
protein expression amounts found in colorectal cancer and 
matched normal tissue, identified that tumour tissue seems to 
be 'primed' and more likely to undergo apoptosis.''^ Similar 
priming effects were previously identified experimentally in 
leukaemias and ovarian cancers by quantifying the effects of 
pro-apoptotic BCL-2 peptides on MOMP and mitochondrial 
depolarisation.^^'^^ Mathematical modelling also allowed to 
confidently identify those colorectal cancer patients that 
responded to 5-FU-based adjuvant treatment within the 
overall patient cohort. Similar translational potential was 
demonstrated by modelling the competence to activate 
effector caspases and therefore to execute apoptosis in stage 
2/3 colorectal cancer patients receiving 5-FU-based adjuvant 
therapy.^"* Tumour tissue in the majority of patients not 
responding to genotoxic chemotherapy appeared incapable 
of undergoing apoptosis execution. Similar results were 
recently also reported in patients suffering from glioblastoma 
multiforme, where results from the mathematical modelling of 
the competence of tumour cells to execute apoptosis 
correlated with progression-free survival times.^^ Protein 
expression patterns investigated mathematically in the 
context of their pathway interplay could therefore serve as 
prognostic or predictive combinatorial systems bio- 
markers,^ '"'^'^"^^ and could provide an alternative when 
classical biomarker identification strategies are unsuccess- 
ful.^^ Importantly, for scenarios in which apoptosis is blocked, 
a systems-level analysis can identify the functional causes for 
cell death resistance and therefore can identify strategies and 
intervention points that re-establish apoptosis competence. 
The feasibility of this concept has been demonstrated in 
recent proof-of-concept studies in which individualised 
predictions on the efficacy of apoptosis-sensitising co- 
treatments, such as BH3-mimetics or lAP antagonists, have 
been generated by modelling platforms. ''^'^'^'^^ Taken 
together, such approaches therefore could provide the 
foundation for the development of a new generation 
of systems-based predictive and prognostic tools that in the 
future may contribute to the personalisation of treatments and 



may allow for patient stratification as part of innovative clinical 
trial design. 



Perspectives for the Integration of Apoptosis Pathway 
l\/lodels into Higher-Level and Multi-Scale Modelling 
Frameworks 

A major challenge for the further development and application 
of apoptosis pathway models is their integration into higher- 
level contexts that reflect organisational structures, processes 
and interactions that go beyond the single-cell level (Figure 4). 
Signalling processes that occur in parallel or subsequent to 
apoptosis signal transduction and that directly or indirectly can 
affect cell death responsiveness may need to be considered. 
For example, the interplay of cellular and multi-cellular 
signalling relevant during cancer immunotherapies may need 
to be considered.- Also, fractions of cells Sying by non- 
apoptotic processes may contribute towards determining 
overall treatment efficacy,^^ and these may also significantly 
alter the tumour microenvironment by promoting immune 
infiltration. Furthermore, activated caspase-3 released from 
dying cells may stimulate growth in surviving tumour cells and 
thereby can drive tumour repopulation,''° a function that could 
be related to the non-apoptotic role of caspase-3 in promoting 
wound healing through prostaglandin E2 signalling.^'' Including 
such paracrine effects and interactions with surrounding 
cells and tissues may therefore be important to assess risks 
for tumour relapse and to more accurately predict short- and 
long-term clinical responses to apoptosis-inducing treat- 
ments. Similarly, the hierarchic organisation of certain tumour 
types, in which cancer stem cells may represent a distinct sub- 
population that exists independent of stochastic cell-to-cell 
differences in gene or protein expression in the overall tumour 
cell population^^'^^ may need to be considered to accurately 
model and simulate tumour growth, invasion and therapy 
response over time. Initial studies in this direction were 
recently conducted and provide data on the behaviour of cell 
populations.^"*'^^ As these population-based models are not 
yet linked with intracellular signalling pathways of apoptosis, 
this may provide an attractive avenue for future research. 

The inclusion of pharmacokinetic (PK) and pharmaco- 
dynamic data (PD) and associated existing PK/PD models 
provides immediate opportunities for higher-level integration. 
Models capable of calculating drug penetration and resulting 
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Figure 4 Systems modelling in the context of biological scale and complexity. 
Biological scale, biological complexity and the associated data availability and data 
quality influence systems modelling strategies. The most insightful single-cell and 
subcellular models of signalling pathways are based on kinetic simulations of 
reaction networks. The complexity of large-scale and population-level data sets are 
mostly tackled by statistical modelling approaches 
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concentration profiles may be of particular importance to 
accurately replicate spatiotemporal profiles of apoptosis- 
inducing stimuli in solid tumours. In this context, tumour 
vascularisation and its modulation by angiogenesis inhibitors 
may also need to be given special attention. First solutions and 
modelling frameworks that allow multi-scale integration have 
now been developed and can be applied in proof-of-concept 
studies.^^ Furthermore, higher-level statistical modelling 
approaches that integrate patient clinical data as well as early 
treatment response data have already been used for persona- 
lising and optimising treatment schedules for chemo- and 
immunotherapies.^^'^^ It is also hoped that ultimately the 
higher-level integration of pathway models will enable complex 
and meaningful mathematical simulations of individual patients 
and patient cohorts before entering clinical trials (Figure 4). 
Such virtual trials could contribute towards assessing treatment 
efficacy a priori and may reduce the currently unsustainably 
high failure rate of phase I I/I 1 1 clinical trials. ''^^^ In addition, the 
need for predictive tools also exists for currently approved 
drugs and treatment regimens, as reliable personalised 
response predictions would reduce the number of unnecessary 
treatments and associated treatment costs. 

Future Challenges and Requirements for (Apoptosis) 
Systems Biology and Systems IVIedicine 

A general challenge both for basic and translational systems 
biological research of cell death processes is the lack of 
quantitative data,^^ in particular at the level of protein 
expression, protein activities, interactions and turnover rates. 
Although several studies made efforts to accurately measure 
the expression amounts of apoptosis-relevant proteins in 
widely used cell line systems, ""^'^^'^^'^^'^^'^^ the determination 
of protein quantities lacks commonly accepted and estab- 
lished methodological standards and may therefore be 
subject to considerable heterogeneity. In the clinical context, 
protein quantity measurements are even less common. 
Routine immunohistochemical analysis in pathology labora- 
tories is limited by low dynamic range and usually not 
calibrated against samples with known protein content. 
Alternative methods that make use of tissue protein extracts, 
such as reverse-phase protein arrays, can provide quantities, 
but lack information on cell-to-cell heterogeneity within the 
tissue context and cannot be used in the routine setting. 
Clinical protein profiling techniques such as multiplex-ELISA 
will require further development to fulfil the needs of systems 
biological and systems medical studies. Although larger 
biobanking initiatives are now starting to widely collect tissue 
for subsequent analysis in research projects, sample storage 
and processing is not uniformly standardised, and the 
available amount of samples are often very limited. This is 
of particular concern regarding the collection of metastatic 
tumour tissue, as most novel therapeutics are initially tested in 
the metastatic setting. Limitations could be alleviated by a 
more comprehensive and earlier integration of clinician 
scientists into the design and coordination of translational 
systems biological research activities, and by the develop- 
ment of novel protein profiling techniques that can be 
integrated into clinical workflows. Furthermore, efforts need 
to be made to disseminate and communicate progress in the 



field of systems biology and systems medicine to medical 
professionals, to expert reviewers of clinical trial protocols and 
to patients and patient advocacy groups. Only in such an 
educated environment the transition from mathematical and 
experimental systems biology towards applied and validated 
novel systems medical approaches can be successful. 
The rapid progress and promising achievements in the field 
of apoptosis systems biology and apoptosis sensitising 
therapeutics may serve as a 'raw model' to facilitate this in 
the coming years. 
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